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SHEAR-STRESS, EDDY- VISCOSITY, AND MIXING-LENGTH DISTRIBUTIONS 
IN HYPERSONIC TURBULENT BOUNDARY LAYERS 

By Dennis M. Bushnell and Dana J. Morris 
Langley Research Center 

SUMMARY 

A method is developed to obtain shear-stress, eddy-viscosity, and mixing-length 
distributions from hypersonic turbulent boundary-layer profile data by using the assump- 
tion of local similarity. Application of the method to published profile data indicates an 
increase in the maximum value of mixing length and eddy viscosity as the Reynolds num- 
ber decreases. The level and trend of this low Reynolds number effect agree with 
similar results previously obtained at low speeds when wall conditions are used in the 
correlating parameter. Examination of the wall-damping effect on the mixing length 
indicates more validity for evaluating the Van Driest damping scale by using wall condi- 
tions than by using local sublayer properties and that the impermeable -wall damping 
constant may be somewhat greater at hypersonic conditions than the usual low-speed 
value of 26. In several instances, areas of disagreement occurred between the expected 
and observed behavior of the data, usually concerning the shape of the published velocity 
profiles. Indications are that the total -temperature data obtained in hypersonic, low 
Reynolds number, cold-wall boundary layers may not be suitably corrected for conditions 
near the wall. 


INTRODUCTION 

Of the several finite -difference and sophisticated integral methods available (refs. 1 
to 11) which compute the development of compressible turbulent boundary layers, all of 
these methods except that in reference 3 utilize either a mixing-length (refs. 1,2,4, 
and 5) or eddy-viscosity (refs. 6 to 11) model for the turbulence-shear term. Most of the 
shear-stress models employed in these calculation procedures were developed from data 
obtained in low-speed flow studies. The application of these low-speed models to com- 
pressible boundary layers was investigated in reference 12, where a set of "universal" 
velocity profiles were derived and used to obtain the turbulent shear-stress distribution 
from the momentum equation. These results indicated that the low-speed mixing-length 
models were directly applicable to adiabatic flows up to a Mach number of 5 and that the 



eddy-viscosity expression of reference 13 could be used by retaining the incompressible 
or kinematic definition of displacement thickness. 

The validity of these low-speed turbulence -shear models for supersonic turbulent 
flows at high Reynolds number was established by comparisons of profile calculations 
with experimental data. (See refs. 1 to 2 and 4 to 11.) Also, work in reference 14 and 
private communication with the author of that reference indicate that the mixing-length 
and "Clauser" (ref. 13) eddy-viscosity models in the outer region are not significantly 
affected by moderate rates of blowing, at least up to a Mach number of 3.5. However, a 
significant increase in the outer-region eddy-viscosity model was found necessary in 
order to successfully compute low Reynolds number supersonic turbulent flows. (See 
ref. 6.) Furthermore, reference 6 points out that as Mach number increases, the low 
Reynolds number type of turbulent flows is more often encountered because of the 
decrease in density and increase in viscosity close to the wall. A similar type of low 
Reynolds number effect (increasing eddy viscosity with decreasing Reynolds number) was 
found for low speeds in references 15 and 16. In addition to this low Reynolds number 
effect, both low and supersonic speed boundary layers show an effect of nonequilibrium 
pressure gradients and high rates of blowing upon the turbulence models. This latter 
effect was pointed out in reference 1, which shows a correlation of the effect of blowing 
and nonequilibrium pressure gradient on the mixing- length model as a function of the 
incompressible or kinematic form factor. The data of reference 17 indicate the existence 
of a nonequilibrium pressure -gradient effect upon the Clauser eddy- viscosity model. For 
supersonic flows, no satisfactory turbulence-shear model exists which includes the com- 
bined effects of low Reynolds number and nonequilibrium pressure gradient. 

Since generally good results were obtained by using essentially incompressible 
turbulent shear-stress models at supersonic conditions, several calculations for hyper- 
sonic conditions have been compared with experimental data. (See refs, 1, 4, 5, 7, 10, 

11, 18, 19, 20, and 21.) These comparisons indicate that the "incompressible" turbulent 
shear-stress models usually give reasonable agreement at hypersonic conditions for sur- 
face phenomena such as heat transfer. However, there are disagreements with the actual 
profile shapes and thicknesses. (See refs. 1, 18, 19, and 21.) 

An alternate method of validating a turbulent shear model (beside comparing com- 
puted and experimental profiles) uses the experimental profile data or correlations thereof 
in the momentum equation and extracts the turbulent shear-stress variation. This shear- 
stress variation in conjunction with the various turbulent-shear-model definitions gives 
the "constants" for a particular model. The shear-stress distribution has been obtained 
by using this method (usually employing velocity-profile correlations) for compressible 
speeds up to a Mach number of 7. (See refs. 12, 14, 22, 23, and 24.) The basic method 
of obtaining mixing-length distributions in this manner was validated for low-speed flows 
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in reference 25. The present work reports on a method to obtain shear-stress, eddy- 
viscosity, and mixing-length distributions directly from hypersonic profile data (Mach 
numbers from 5 to 20) by using the assumption of local similarity and compares results 
by this method with values and trends of turbulent shear-stress models successfully 
employed at low and supersonic speeds. 

SYMBOLS 


Values are given in both SI and U.S. Customary Units. The measurements and 
calculations were made in U.S. Customary Units. 


A 

A + 

Cf 

K 


damping scale, 


A + M 

P u T 


Van Driest damping constant (see eq. (9)) 


skin-friction coefficient, 


w 


1 2 
fPeV 


"Prandtl wall slope," slope of mixing-length variation with y in region 
~ ~ 0.1 and passing through the origin 


l 


mixing length 


M 

N 


Mach number 


exponent of power-law velocity profile, 


u 


\1/N 


n 


integer 


P 

R 

R 

R; 


si 


pressure 


Reynolds number 


Reynolds number based on sublayer thickness, 


Reynolds number based on momentum thickness, 


Pw u e y s i 
Pw 
P e u e 0 


Mt 


radius 


temperature 
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u,v 

longitudinal and normal velocity component, respectively 

U T 

friction velocity, 

u' 

velocity fluctuation 

x,y 

longitudinal and normal Cartesian coordinate, respectively 

1°M W 

sublayer thickness, 

u r,vrw 

y s i 

y+ - yl w 

|w 

p w 


5 

density thickness of boundary layer, value of y where — = 0.995 

Pe 

5 + 

u r w p w 6 

boundary -layer thickness in law of the wall coordinates, — — 

P w 

6* 

displacement thickness 

R* 

5 i 

incompressible or "kinematic" displacement thickness 

e 

kinematic eddy viscosity 

e 

momentum thickness 

M 

molecular viscosity 

p 

density 

T 

Ojl 

total shear stress, (jit + pi) — , unless otherwise noted 

Subscripts: 

« 

e 

local external to boundary layer 

max 

average maximum value 

o 

stagnation condition 
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t turbulent 

w wall 

ANALYTICAL APPROACH 


Basic Equation 

The equations for the conservation of mass and momentum in axisymmetric or two- 
dimensional compressible turbulent boundary-layer flows (see ref. 8, for example) are 

Continuity: 

9(pur) 3(pvr) _ 0 ' (1) 

9x 9y v ’ 

Momentum: 



" r + ( rT ) 
ax ay v 


(2) 


where, for internal or nozzle-wall flows (the usual case of available hypersonic data, 
ref. 26) r « r w - y. The shear stress includes the turbulent component and is defined 
herein as r = (p + pe) Note that the boundary-layer thickness can be the order of 
the transverse radius of curvature. 

An expression for pvr is obtained by integrating equation (1) once with respect 
to y. By assuming an impermeable wall, the following equation is obtained: 



The usual assumption is next made that u/u e , p/p e , and p/p e are all functions of only 
y/5 in the neighborhood of the profile under consideration. (See ref. 22.) The 5 value 
used herein is the density thickness, which is perhaps more meaningful at hypersonic 
conditions than a velocity thickness. The velocity profiles of reference 21 at M e = 20 
are shown for various stations along the nozzle in figure 1 and indicate that for these 
data the assumption of a "local similarity” in y/6 coordinates is acceptable, even at 
M e = 20. Obviously, in severe nonequilibrium flows, the assumption of y/6 similarity 
would not be valid, even locally. Therefore, the data examined will be restricted to 
hypersonic flows where the local pressure gradient is small or zero (although the nozzle - 
wall boundary layers have a favorable pressure-gradient history). It should be noted 
that the nonsimilar terms could still be important in the near-wall region, as the sub- 
layer thins in comparison with 6 with increasing Reynolds number. 
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By using this assumption ^u/u e , p/p e , p/p e are a U functions of y/5), inserting 
equation (3) into equation (2), and integrating equation (2) out from the wall,, the following 
equation is obtained ^with v w = 0 and noting at high Mach number that pressure is a 


function of y/6 even for = 0, ref. 21^: 


t 'i.iw5 + 


d (Pe u e 2 ) 


+ d5 + 5 dr w 


p e u e 2 ' dx dx r w dx 


"r^yjL/ifd 

J 0 6P e lu e ] 


,2 d 


(Pe u e 2 ) 


r w P e V 


25_d6 
dx r w dx 


f y/5 Ji.idfl') 

^0 u e 


u 6 d (p e u e 

Li- 

d6 

. u_ 

6 

dr w 

u e P e u e dx 

u e 

dx 

u e 

r w 

dx 


r y/5 y p u d /y\ 

u 6 2 d ( p e u e) U 5 2 d6 

J 0 6 P e u e \ 5 ) 

u e r w p e u e dx u e r w ^ 


6 d P t 


Pe Ue 


2 dx 


’f y/ V d fL 

_ S2_dp^ 

f y/6 y p d /y\ 

J 0 p 6 W 

dx 

J 0 6 P e \EJ 


P u ‘ 
H e e 


p e d6 P y , p e d5 

r 


p e u e 2 ^ p e 6 P e u e 2 dx 


r y/5 jl i[i 

J 0 p e u 



5 p e d6| 

ty\ 2 P 

8 p e dS 

( 

■ y/6 P 2y | 

f l) 

is 

r w p e u e 2 dx 

Vdj P e 

r w p e u e 2 dx 

J 


Kd) 


(4) 


Except for reference 21, the data considered herein were reduced by assuming ^ = 1 
for all y; therefore, the sum of the last four terms in equation (4) would not appear. As 
stated previously, this equation includes transverse curvature effects for internal or 
nozzle-wall flows as well as pressure -gradient effects. Two-dimensional flow is recov- 

dr w / 9 

ered by letting r w — °° and — 0. Equation (4) constitutes a relation for T /P e u e 
once the profiles of u/u e , p/p e , and p/p e as functions of y/5 and the various 
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x-derivatives are known. In the present work, these profiles and derivatives are taken 
from experimental data. Therefore, accurate profiles and values of p Q and u e at 
several longitudinal stations are necessary for determination of the x-derivatives. 

Because of the form of the basic equations (eqs. (1) and (2)), equation (4) cannot be 
used for flows with significant longitudinal curvature. The pertinent equations for this 
case are available in reference 27, and these would have to be integrated in a manner 
similar to that used herein to handle flows with significant longitudinal curvature |where 

JL > — L- for the turbulence models, ref. 28). With the shear-stress distribution calcu- 
r w 300 J 

lated from equation (4) , shear stress is obtained from the relation 

a 8 ( u / u e) (5 

P e u e 2 P e u e 2 P e U e 6 9 ( y / 6 ) 


The definitions of the various turbulent shear-stress models are used with 
to determine the eddy-viscosity and mixing-length distributions; that is, 


r t/ 


P e u e 


and 



( 6 ) 


(7) 


which were obtained from the definitions 


r t = pe 


9u 

9y 


and 


r t = P l 2 


9u 

9y 


8u 

9y 


The derivatives of the experimental velocity profile 


8 ( u / u e) 

9(y/6) 


are needed to obtain e 


and Z variations with y/6. Also, when the near-wall data are sufficiently detailed to 
allow an examination of wall-damping effects, the expression from reference 29 (Van 
Driest damping function) will be employed: 


where 


f=f(i-e-y/ A ) 


( 8 ) 

(9) 
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and A + is the usual "damping constant" (a + = 26). At the present time, no clear choice 
exists as to whether the p, p, and r in equation (9) should be evaluated at the wall or 
locally as a function of y. This question is investigated in the present paper. 

In the present method, equation (4) is used to obtain shear-stress, eddy-viscosity, 
and mixing-length distributions from experimental data with the following restrictions: 

(1) Impermeable wall 

(2) Small or zero streamwise pressure gradient 

(3) Accurate profile data at several longitudinal stations 

(4) Small longitudinal curvature 

(5) Mach numbers greater than approximately 5 (hypersonic) 

Unfortunately, the amount and, to a certain extent, the accuracy of the published data 
decrease significantly from low-speed to supersonic and then hypersonic conditions. The 
data currently available which fulfill the restrictions noted are basic hypersonic turbulent 
boundary -layer information, and such data are not plentiful; however, the data presented 
in references 21 and 30 to 36 will be used herein to exercise the present method, and the 
results should give some indication of what shear-stress models might be used in hyper- 
sonic flow. 


Method of Solution 

Profiles of u/u e and p/P e as functions of y/6 were plotted for the data from 
references 21 and 30 to 36. Curves faired through each data point (except as noted in 
the discussion of the individual cases) were used to obtain detailed p/p e and u/u e 
distributions for the integrals in equation (4). In order to achieve a high degree of accu- 
racy, 200 points were read from each curve and used as input to a computer program in 
tabular form. The trapezoid rule was used to evaluate the integrals in equation (4) with 
the program interpolating four points between each pair of points in the input table; this 
interpolation gave a 1000-point integration for the entire boundary layer. This integra- 
tion allowed the determination of r/p^u^ as a function of y/6. The necessary local 
streamwise derivatives were obtained graphically from fairings of the quoted variations 
o f 6, r w , p e u e 2 , p e u e , and p 0 . 

The derivatives of u/u e with respect to y/6, necessary to obtain eddy-viscosity 
and mixing-length values from the shear-stress profile, were computed by using a least- 
squares method with an overlapping series of curve fits for six consecutive points in the 
input u/u e table. Each curve fit was differentiated, and the derivatives of the third 
and fourth points in the curve fit were computed from this differentiation. The third and 
fourth points used in the first curve fit were then taken to be the first two points in the 
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next curve fit and so on. The derivatives obtained by this method were plotted and 
faired. The resultant derivatives are believed accurate to within approximately 15 per- 
cent (assuming the experimental data are exactly correct, which may not be true in all 
instances). The velocity derivatives in the outer region were successfully checked by 
comparison with results from an "N -power", approach; that is, a velocity distribution 

— = was assumed and differentiated. By using the experimental N value, a 

u e \o/ 

local derivative in the outer region was obtained to check the results from the curve-fit 
procedure. 


SHEAR-STRESS, EDDY- VISCOSITY, AND MIXING-LENGTH DISTRIBUTIONS 

The results of evaluating equation (4) along with the e/u e 5* and l/5 distribu- 
tions for typical profiles from references 21 and 30 to 36 will now be presented. The 
various x or longitudinal derivatives and other input and output quantities used in the 
solution are given in table I. The skin-friction values used were obtained wherever pos- 
sible from values quoted in the reference from which the data were taken. Exceptions to 
this procedure are noted in the discussion of the individual profiles considered. It should 
be noted that because of the incomplete nature of much of the published hypersonic data, 
the results presented herein of applying the present method should be considered as indi- 
cating trends rather than true absolute levels. 

Mixing-length and eddy-viscosity distributions from reference 12 will be shown for 
comparison with the present computed hypersonic results. The curves shown from ref- 
erence 12 are for supersonic flows over a fairly wide range of Reynolds numbers and are 
typical of variations currently employed in prediction methods. 


Axisymmetric-Nozzle-Wall Data at a Mach Number of 20.8 


Profile data are available from reference 21 for several stations along an axisym- 
metric nozzle wall at M e ~ 20 in helium. The wall-to-total temperature ratio was 
near 1.0, and the data were reduced by assuming T 0 was constant. The large density 


ratio across the boundary layer — = 140 provided a severe test case of compressibility 

T - T /„ \2 

or high Mach number effects. A re-reduction of the data by using 0 w 


To,e - T 


w 


(k) • with 


the lowest possible T w /T 0 being 0.96, indicated negligible change on the velocity profile 
for ^ > 0.1. Results are shown in figure 2 for both total-temperature distributions. 


The profile chosen was from station 108 with a local Mach number of 20.8 and a 
small local favorable pressure gradient. (See table I.) The profile is given in table 2 
of reference 21. The p/p e profile used in equation (4) was the same one used to reduce 
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the profile data in reference 21. The Rp was 6.90 x 10^, which at supersonic speeds 
is usually large enough to insure a fairly well-developed turbulent flow. (See ref. 37.) 

The shear-stress results are shown in figure 2(a). As stated previously, the 5 used 
in the present work is the pitot or density thickness. The nominal velocity thickness is 
less (—= 0.995 corresponds to ^ ~ 0.5 V The computed shear stress goes to the lami- 

\ u e 6 / v 

nar value when — = 0.9983 and *=■ = 0.695. The fact that r f — 0 at ~ < 1 may be 

u e 6 1 5 /o/, \ 

/ a(u/u e ) 

reasonable in this case because of the extreme fullness of the velocity profile — - 7 — 0 

\ \9(y/6) 

when ^ < 1). Also, the hot-wire data of reference 21 for total-temperature fluctuations 
may be interpreted in the outer region as twice the velocity-fluctuation level because of 
the high Mach number, and these results indicate a u'/u e level less than 0.5 percent 
for ^ > 0.7 (low vorticity level). 

At high Mach numbers, the energy equation reduces to u thus, for a hot- 

or near-adiabatic-wall turbulent boundary layer at high Mach number, the u/u e profile 
is extremely full because T 0 is nearly constant. However, at the same Mach number 
and a cold-wall case, where T c may decrease in the outer high Mach number portion 
of the boundary layer, the velocity profile would be considerably less full. (See refs. 37 
and 38.) Therefore, at hypersonic speeds, another mechanism besides the turbulence 
field is available which can have a first-order effect upon the mean velocity profile. 

This discussion implies that at hypersonic speeds correlations of pressure gradient and 
wall-blowing effects which are based on velocity- profile shape alone (such as those given 
in reference 1) may not be reliable because the correlations are based upon low-speed 
data where the turbulence is the dominating factor on the profile shape. 

The eddy-viscosity and mixing-length distributions are shown in figure 2(b) and 2(c) , 
respectively. In the outer region, the e/u e 6* and 1/6 values are significantly larger 
than the usual adiabatic equilibrium value of 0.018 and 0.09, respectively. (See ref. 12.) 

A later portion of the present paper shows that these higher values may be expected and 
are a logical extension of the low Reynolds number effect noted in references 15 and 16 
at low speeds. As indicated in the figure, the wall-region < 0.3^ results in figure 2(c) 
indicate a Prandtl wall slope of approximately 0.65 rather than the more usual value of 0.4 
(ref. 12). This increased value of the slope again is shown later to be an expected result 
of the low Reynolds number trend noted in reference 39. Further calculations using the 
method of reference 1 indicate that the larger Prandtl slope at the wall constitutes the 
major reason for the disagreement shown in reference 21 between these profile data and 
predictions using the method of reference 1, which originally employed K = 0.4 . 

The mixing-length distribution shown in figure 2(c) indicates a region & < 0.1^ where 
the Prandtl wall variation ( l = Ky) is damped by the wall. This damping effect will be con- 
sidered in a later section. 
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Turbulent Boundary-Layer Development Inside a Circular Duct 

at a Mach Number of 6 

Profile data for the development of a turbulent boundary layer inside a constant- 
area circular duct at a Mach number of 6 are presented in reference 30. The duct was 
connected to the exit of a contoured axisymmetric nozzle, and the measurements were 
made during the subsequent development of the nozzle-wall boundary layer. Because of 
the 5* increase in the duct, there is a small positive pressure gradient, whose strength 
is a function of Reynolds number. Two profiles which are typical of the low and high 
Reynolds number range of the tests were analyzed. The first profile examined was 

obtained at station 124, p Q = 0.45 MN/m 2 (65.3 psia), — ^ = 0.73, and M e = 5.76. The 

T o 

profile is given in table in(a) of reference 30. These results are typical of the low 
Reynolds number range of these data (Rq = 8 x 10^). 

The shear-stress distribution is shown in figure 3(a). Here the shear does not fall 
monotonically from the wall value but peaks at ^ ~ 0.2 before falling off toward the outer 
edge. Typically, this peaking away from the wall occurs for boundary layers in the pres- 
ence of adverse pressure gradients or wall blowing. (See ref. 40.) Note that the velocity 
profile is much less full than for the Mach 20 helium case (shown in fig. 1) and exhibits a 
more usual shape (N = 7.7). The eddy-viscosity and mixing-length variations shown in 
figure 3(b) and 3(c), respectively, are fairly similar to those observed at M e = 5 and 
high Reynolds number. (See ref. 12.) Therefore, since the Rp for this profile 
(8.00 x 10 2 ) is fairly close to the value for the Mach 20.8 profile (6.90 x 10 2 ) and the l/5 
and e/u e 6* profiles are quite different (compare figs. 2(b) and 2(c) with figs. 3(b) and 
3(c)), it appears that the magnitude of Rg is no longer a valid representation of the 
"relative turbulent state” of a boundary layer at hypersonic speeds but the turbulent state 
is a function of Rg, M e , and T w /T 0 . For the Mach 20.8 profile, the wall-to-edge den- 
sity ratio = 0.007j is much less than that for the Mach 6 profile = 0.18^. This 

difference in the P w /P e ratio represents a probable explanation for the discrepancy. 
Also, the u T value is less at Mach 20.8 because of the much lower Cf value. There- 
fore, at hypersonic speeds a Reynolds number based on wall rather than edge conditions 
probably better characterizes the structure of the boundary layer. It is shown later in 
the present paper that the use of wall parameters allows a reasonable correlation of the 
available low Reynolds number eddy-viscosity and mixing-length data to be made from 
low speeds up to a Mach number of 20, at least for near-adiabatic flows. 

The other profile examined from reference 30 was at station 172 and 

T 

p Q = 2.17 MN/m 2 (315 psia), —■ = 0.7, Rg = 47 x 10^, and M e = 5.75. The profile 

was taken from table IV(c) of reference 30. The shear-stress profiles are shown in 
figure 4(a). For this profile, the shear did not reach zero at the edge but deviated at 
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y 

^ « 0,75. The reason for this deviation is not known, but probable explanations are 
(1) y/5 similarity assumption erroneous or, perhaps more probable, (2) incorrect val- 
ues of dp/d x and dS/dx. The profile measuring stations are fairly far apart, which 
makes the evaluation of accurate local x-derivatives difficult. A curve has been shown 
drawn to zero at the edge, and results will be shown for both computed and assumed 
outer-shear distributions. Again, the peaking away from the wall is apparent. A com- 
parison of figures 3(a) and 4(a) indicates that the near-wall region (where significant lam- 
inar shear occurs) constitutes an appreciably smaller fraction of y/ 6 in figure 4(a) as 
would be expected because of the higher Reynolds number. The eddy-viscosity and 
mixing-length distributions are shown in figures 4(b) and (c). At this higher Reynolds 
number, the values are still in reasonable agreement (although somewhat higher than 
those at lower Mach numbers) with the results from reference 12 for M e i 5. This 
agreement corresponds to the finding in references 15, 16, and 39 that at high Reynolds 
number the 1 / 5 , e/u e 6*, and K values are nearly invariant with Reynolds number (the 
cornerstone of most finite-difference calculation methods). 

Conical-Nozzle-Wall Data at a Mach Number of 10.65 

Profile data for the development of a turbulent boundary layer on a conical nozzle 
wall (favorable pressure gradient) are available in reference 31. The nozzle had a total 
conical expansion angle of 7.5° and was attached to a hypersonic gun tunnel facility. In 
the data region, the nominal test-section diameter was 0.102 meter (4 inches). Over a 
distance of 0.204 meter (8 inches), the Mach number increased from 10.2 to 11.6. The 
ratio of wall to total temperature was low (0.26); therefore, especially accurate total- 
temperature profiles are necessary in order to correctly reduce the pitot data to velocity 
profiles. (The static pressure was assumed constant.) The data analyzed herein cor- 
respond to x = 0.665 meter (26.2 inches), M e = 10.65, and p Q = 30.7 MN/m^ (4450 psia). 
The profiles were not tabulated and hence had to be obtained from figures 5.10c and 5.11b 
of reference 31. The Rg was 12.8 x 10^. This profile is typical of the high Mach num- 
ber, high Reynolds number (comparatively speaking) portion of the data in reference 31. 

The shear-stress results are shown in figure 5(a). The lower set of curves cor- 
respond to using a Cf obtained from table 5.1 of reference 31. A value of 5.2 X 10"^ 
was taken as a mean of the "measured" results, which varied from 5.49 X 10"^ to 
4.89 X 10" . This value gives a shear stress which approaches zero at ^ = 0.5, clearly 
not correct for a velocity profile which has an N power of 6.2 (ref. 31). Solutions for 
shear stress obtained by varying the input parameters over a reasonable range of uncer- 
tainty yielded shear-stress distribution not very different from that shown. The only 
perturbation on the inputs which would give a reasonable shear-stress distribution was to 
increase the Cf value. Results are indicated by the upper set of curves for an arbitrary 
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increase of 100 percent in Cf. These results do not imply that the actual Cf for the 
profile was as high as 1.04 x 10“^, but rather that something is wrong with either the 
present approach and/or the profile data. The validity of the local similarity assumption 
^u/u e and p/p e are functions of y/s) can be questioned; however, the large effect 
shown in figure 5(a) indicates that there may be some problems with the data, particularly 
with the total -temperature profiles. 

Reference 1 indicates that to compute this profile in terms of a total -temperature— 
velocity relation, the total turbulent Prandtl number (see ref. 1) would have to be increased 
to 2 in the outer region of the boundary layer. (See fig. 6 of ref. 1.) This value is con- 
siderably above the usual values of total turbulent Prandtl number (near 1); therefore, the 
total-temperature data for < 0.8 ^near the wall, ^ < 0.2^ may be too low. This 

assumption may be reasonable in light of the short running time of the facility (10 to 
20 msec), the existence of a low Reynolds number region near the wall where the T 0 
probe corrections are presumably different from those obtained in the free stream (see 
ref. 38), and the fact that numerous runs with slightly varying conditions were necessary 
to obtain a single T 0 profile. Consequently, the extreme difficulty of obtaining accurate 
T 0 measurements near the wall in such a facility suggests that for ^ < 0.2 the T 0 
measurements may be low. This inaccuracy could explain part of the problem with the 
shear stress approaching the laminar value at ^ = 0.5. If the T 0 values near the wall 
were higher, the u/u e profile would necessarily be fuller (higher N value); therefore, 
the integrals in equation (4) would be affected. Also, the indicated Cf would probably 
increase. The net result of higher T 0 values would be some increase in Cf and a 
change in slope of the curve of u e ^ as a function of y/6. 

Eddy-viscosity and mixing-length distributions are shown in figure 5(b) and (c), 
respectively, for the two r/p u e ^ variations of figure 5(a). Using the upper shear- 
stress curve from figure 5(a) results in more conventional values of l/5 and e/u e 5f 

for this fairly high Reynolds number boundary layer. However, if the T 0 profile was 

3(u/u e ) 

fuller, thus increasing the N power on the velocity profile, the — / values would be 

9 (y/ o) 

lower in the inner region of the boundary layer; therefore, the increase of l with y 
would be greater for ^ < 0.2. This discussion would explain the low Prandtl slope for 
this case (0.26 versus usual values of 0.4 or greater). 

Axisymnietric -Nozzle -Wall Data at a Mach Number of 12 

Profile data for the development of a turbulent boundary layer on a contoured axi- 
symmetric nozzle wall with a mild favorable pressure gradient are given in reference 32; 
tabulations of the profiles were supplied by the author of that reference. The profile 
examined herein was taken at station -0.026 meter (-1.023 inches) (measured from the 
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nozzle exit) on a contoured steel nozzle with p Q = 6.9 MN/m 3 (1000 psia). The local 

T 

Mach number was 12.0 and —■ = 0.295. 

T o 

The scale of the facility is similar to that used in reference 31, with a nominal 
test-section diameter of 0.102 meter (4 inches). However, the Mach number is higher 
(12 versus 10.6) and the total pressure lower (6.9 MN/m 3 (1000 psia) versus 30.7 MN/m 3 
(4450 psia)); therefore, the Rg value of 1.12 X 10 3 is considerably below the value for 
the profile examined from reference 31. In fact, an Rg. value of 1 x 10 3 is low even 
for a supersonic boundary layer. (See ref. 37.) It should be noted that this profile con- 
stitutes some of the highest Reynolds number data taken in the investigation of refer- 
ence 32. The possibility that this represents an "overdamped" boundary layer will be 
discussed in a later section. 

The shear-stress profiles are shown in figure 6(a). Again, as in figure 5(a), the 
Cf value had to be increased by a factor of nearly 2 in order for the shear stress to 
approach zero in the outer region of the boundary layer. However, such a higher Cf 
value is perhaps defensible from the velocity-profile data shown in figure 6(b). These 
profile data raise an interesting point, one brought out in reference 33, as well as else- 
where. The data points near the wall in figure 6(b) exhibit a deviation from the general 
trend of the data. This deviation might be attributed to probe-interaction effects with the 
wall. However, another possible cause for at least some of this deviation is indicated in 
figure 16 of reference 33. In the very-near-wall region, especially for this low Reynolds 

/y,i \ 

number, turbulent boundary layer with a thick sublayer I — — = 25 percent j, the laminar 
shear stress dominates a significant portion of the inner profile (out to y/6 from 0.05 
to 0.2 from fig. 6(a) depending upon the Cf value). Therefore, the shear stress, 
which is nearly constant, is given by ju. -|^. 

Because of the heat transfer to the wall, the static temperature (and conse- 
quently fj.) increases to a peak^at | » 0.12, ref. 32*j before decreasing to the edge value. 
The requirements of t being constant and \i peaking imply that 9u/9y must first 
decrease and then increase again as the distance from the wall increases. Therefore, 
the velocity profile should exhibit a "kink" as illustrated in figure 16 of reference 33. 

The direction of such a kink is consistent with the deviation of the first data point in fig- 
ure 6(b) from the rest of the profile. The consequence of this discussion is that the avail- 
able data near the wall may not be sufficient to define a true Cf value; therefore, a 
higher Cf (such as used in fig. 6(a)) may be reasonable for this boundary layer. Fur- 
ther research is obviously needed to clarify the phenomena in the wall region of a low 
Reynolds number, hypersonic turbulent boundary layer, especially under cold-wall 
conditions. 
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The eddy -viscosity and mixing-length results are shown in figure 6(c) and 6(d), 
respectively. The values of both l / 6 and e/u e 6^ are rather low; this result will be 
discussed more fully in a later section. 


Data Obtained in a Two-Dimensional Boundary-Layer Channel 
at a Mach Number of 5 


Profile data are available from reference 33 for the development of a turbulent 
boundary layer at a Mach number of 5 in a two-dimensional boundary- layer "channel" 
with essentially zero pressure gradient locally. The channel wall is a continuation of 
a Mach 5 two-dimensional nozzle; therefore, the data were essentially from a two- 
dimensional nozzle-wall boundary layer, which is allowed to "relax" in a zero pressure 
gradient. Two profiles corresponding to runs 6211 and 12091 are examined. (See 
table l(r) and 1(d), respectively, from ref. 33.) The local Mach number is 4.7 and 



“ q 

The first profile examined (run 6211) corresponds to an Rg = 57 x 10° and, there- 
fore, should constitute a high Reynolds number flow. The shear-stress, eddy-viscosity, 
and mixing-length distributions are shown in figure 7. The shear-stress distribution 
looks fairly reasonable; however, the eddy- viscosity and mixing-length distributions show 
a "trough" region for 0.1 < ^ < 0.6 approximately. In this region, both eddy-viscosity 
and mixing-length profiles are not as full as one would expect, especially for this Mach 
number and Reynolds number. (Note comparison with M e ^ 5 results from ref. 12.) 


A possible reason for the disagreement with the reference 12 results is readily 
available from the velocity profile for this case (shown plotted in log-log form in 
fig. 5c(2) of ref. 33). The N-power quoted in reference 33 for this profile is 8.3. Fig- 
ure 5c(2) of reference 33 indicates that this N-power is valid for y/6 between 0.02 
and 0.1 and also between 0.6 and 0.98. However, for y/5 from 0.1 to 0.6, the indicated 

9 ( u / u e) 


N-power is significantly lower (N ~ 5). Therefore, there is an increase in 


9(y/s) 


m 


the midregion of this boundary layer which is not normally expected and which causes 
"less full" eddy-viscosity and mixing-length profiles in this region. 


This apparent abnormal behavior of the variation of u with y may be due to 
(1) residual nonequilibrium effects of the nozzle expansion upon the boundary layer (see 
ref. 26), (2) nozzle sidewall effects, or (3) inaccurate total-temperature profiles. How- 
ever, how these could cause this anomalous behavior in the velocity-profile data is not 
clear. The eddy-viscosity and mixing-length values which will be used in a later section 
of the paper will correspond to the more "normal” sections of the profile < 0.1 and 

i>°4 
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Results from the second profile (run 12091) are shown in figure 8. The Rg for 
this profile is 5 x 10 , which is still comparatively high but an order of magnitude less 
than for the previous profile. Again, the shear-stress results are reasonable and the 
1 / 5 and e/u e 6^ distributions have the same odd-looking trough as the higher Reynolds 
number data. This region has roughly the extent 0.1 < ^ < 0,65. In figure 5b(l) of ref- 
erence 33 (a log-log plot of velocity as a function of y/6 of this profile), there is again 

y 

a corresponding "low N" region for 0.1 < ^ < 0.65. Therefore, at an order of magni- 
tude lower Reynolds number, the same anomalous phenomenon has occurred at roughly 
the same location. 


P /Pe 

used: 


To check the validity of the "local similarity" assumption for this case (u/u e and 
are functions of y/s), a simple form of the integral momentum equation can be 


£f 

2 



0 d6 
5 dx 



( 10 ) 


For a flow in local similarity, the ratio 0/6 should be invariant with x; there- 
fore, the second term in equation (10) represents a nonsimilar correction. For the 

C f 4 

first profile considered (fig. 7), the terms in equation (10) are — - = 3.6 X 10“’ and 

Cl 

3.86 x 10“^j therefore, the second or nonsimilar term in equation (10) constitutes 

only about a 7-percent correction, and this flow is indeed near a local similarity, at least 
for the profile taken as a whole. 


Measurements on a Hollow Cylinder at a Mach Number of 6 

References 34 and 35 provide profile data for a boundary layer grown on a hollow 
cylinder at a Mach number of 6. Reference 34 presents an adiabatic-wall case, whereas 
reference 35 shows a cooled-wall condition on the same configuration in the same facility. 
The pressure gradient remains close to zero throughout the entire history of the boundary- 
layer development; therefore, these data should be "clean”; that is, there are no questions 
concerning possible effects of upstream pressure-gradient history, as was the case for 
the nozzle-wall results considered previously (figs. 2 to 8), However, as was shown in 
reference 37, the presence of a transition region provides a perturbation in the profile 
shape which requires a hundred or so boundary-layer thicknesses for relaxation. The 

profile examined from reference 34 is tabulated in table 1(g) of that reference and cor- 

T 

responds to station 0.15 meter (6 inches), M e = 6.02, — 0.895, and R$ = 1.4 X 10^. 

T 0 

The shear-stress distribution is shown in figure 9(a), and the eddy-viscosity and 
mixing-length distributions in figure 9(b) and 9(c), respectively. The eddy- viscosity 
distribution peaks somewhat inboard, as indicated by comparison with the reference 12 
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results. The mixing-length distribution exhibits a high K value near the wall but a 

modest level in the outer region. Generally, a high K value is accompanied by a rather 

high level of (4) .A possible reason for this behavior was reported in reference 37, 
Wmax 

which shows that these profiles have a rather high N -power (9) for this Reynolds number. 
The profile evidently has not relaxed completely from the "N -power overshoot" (extreme 
profile fullness) that seems to dominate regions immediately downstream of transition at 
hypersonic speeds. (See ref. 37.) 

T w 

The second profile examined was from reference 35, where — — = 0.45, 

A T ° 

R q = 1.28 x 10 4 , and M e = 5.98. The profile is from station 4 and is shown in table 1(d) 
of reference 35. The shear-stress, eddy-viscosity, and mixing-length distributions for 
this profile are shown in figure 10. The eddy-viscosity (fig. 10(b)) and mixing-length 
(fig. 10(c)) distributions are somewhat lower than those of reference 12. 

Measurements on a Flat Plate at a Mach Number of 6.5 

Reference 36 provides data for a boundary layer developed on a sharp flat plate at 

T w o 

M e = 6.5. The profile examined herein corresponded to = 0.48, Rg = 2.53 x 10 , 

T 0 

and zero blowing with the test designated as run 3 and the profile given in table 5 of 
reference 36. 

The shear-stress, eddy-viscosity, and mixing-length distributions are shown in 
figure 11. The shear-stress profile was computed by using a value of Cf/2 of 
9.5 x 10” 4 for this case, as used in reference 24 instead of 7.23 X 10 -4 , as used in ref- 
erence 36. The eddy- viscosity values are somewhat above the results from reference 12. 
As discussed in connection with previous profiles, this result is perhaps reasonable 
because of the low Reynolds number involved. The wall-damping region seems to extend 
out to ^ ~ 0.15. Therefore, the conventional region of the Prandtl wall slope is rather 
abbreviated, and the mixing levels off at a reasonable value. (See fig. 11(c).) 

COMPARISON OF HYPERSONIC AND LOW-SPEED EDDY-VISCOSITY 
AND MIXING-LENGTH LEVELS 

As stated previously in the present report, a definite low Reynolds number effect 
was noted in references 15 and 16 on the / — and (-] levels in low-speed 

V u e 6 i/ max W max 

flows. A similar effect was observed for K (increase with decreasing R) in refer- 
ence 39 for low-speed flows. All of these investigations (refs. 15, 16, and 39) utilize 
R q as the correlating parameter for low Reynolds number effects at low speeds. At 
hypersonic conditions, one expects a similar low Reynolds number effect. However, 
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comments in the present paper in connection with the first two cases considered indicate 
that at hypersonic conditions (where generally p w « p g and /i w » p e ) Rg, which is 
based on edge conditions, would not be expected to correlate low Reynolds number effects. 
This disagreement between low-speed and high-speed data when Rg is used as the basis 
for comparison is indicated in the inserted figure from reference 39 in figure 12(a). 
Reference 41 obtained good results in correlating high Mach number (M e = 20 to 43) pro- 
files by using wall conditions and a wall-friction velocity. Therefore, in an attempt to 

correlate low Reynolds number effects at both low and hypersonic speeds, a Reynolds 

p u 6 

number based on wall conditions is used. The form chosen was 5 + = - J ■. Another 

Pw 

kinematic scale factor, such as 6?, could perhaps have been used; however, the use of 
5 + allows the ready identification of extremely low Reynolds number flows which may be 
"overdamped," since the wall-damping function (eqs. (8) and (9)) contains a damping con- 
stant A + which is defined in terms of law of the wall coordinates as 
(for impermeable-wall condition). 

Generally, when y + » A + , the profile data are in agreement with some form of the 
logarithmic (or "wall") law; the logarithmic region is considered the fully turbulent por- 
tion of the profile. For y + A + the wall damping affects the profile directly. As 5 + 
approaches A + , the wall damping will affect more and more of the profile and eventually 
reduces the turbulence in what is usually considered the outer region of the boundary 
layer 


A ^ -LZ = 26 = A + 
Pw 


0.2^. (See also reverse transition discussion in ref. 42.) 
The low-speed results of references 15, 16, and 39 for K, ^ 


max 


and 


u e 6 i 


max 

have been plotted as a function of S + (the tabulation of ref. 43 was used to convert from 
Rg to 5 + for the ref. 16 and 39 results) and are shown in figure 12. Also shown are 
the results presented herein (figs. 2 to 11) and those from references 12 and 13. Where 
anomalies have been found in the data (these are discussed in connection with the indi- 
vidual cases), the results have not been shown. The values shown from the results pre- 
sented herein are tabulated in table I. 


The variation of Prandtl wall slope K with 8 + is shown in figure 12(a). The K 
values are only quoted for profiles where the damping does not extend beyond Z « 0.15; 
that is, where a definite Prandtl wall region is present. If damping persists out farther 
than ^ « 0.15, a true law of the wall region probably does not occur and the profile goes 
almost directly from a sublayer to a wake region. Because of possible inaccuracies in 
the present analysis procedure (such as the similarity assumption) and the incomplete- 
ness of some of the data, considerable scatter exists. (This scatter is shown in all parts 
of figs. 12 and 13.) However, the present hypersonic values (especially the data at near- 
adiabatic -wall conditions) do tend to exhibit the same trend as the incompressible values 
of reference 39; that is, the values increase as Reynolds number decreases. Interestingly, 



these K values at low 6 + are appreciably above the usual accepted values of 0.4 
to 0.44. These higher K values have an appreciable effect upon profile shape and pre- 
dicted Cf when employed in a finite-difference procedure. 

The outer mixing-length levels are shown in figure 12(b), and again the present 

hypersonic results are in fair agreement with the level and trends of the low-speed data. 

However, at very low values of 6 + (66), where 5 + => A + , the (!•) from refer- 

V o/max 

ence 32 has reversed the trend of the values at higher 8 + . This reversal in trend may 
be due to excessive wall damping, since 5 + /a + has a value of only 2.5. 

Obviously, with 6 + sufficiently reduced the flow will revert to a laminar condition. 

(This relaminarization occurred in ref. 44.) Therefore, somewhere the upward trend of 

the low-speed data has to reverse, and the values decrease toward zero. The 

' o/max 

actual value of 6 + where this reversion occurs is probably a function of the flow history, 
local conditions, and wall -to-total -temperature ratio and, therefore cannot be correlated 
on a plot such as figure 12(b). 

It should be noted that the bulk of the published hypersonic turbulent profile data 
(such as refs. 21, 30, 31, 32, and 33) was obtained on nozzle walls (see also ref. 26) pri- 
marily because of the difficulty of obtaining turbulent boundary layers on basic configura- 
tions in hypersonic flows. This difficulty can be attributed to the trend toward high tran- 
sition Reynolds number at high Mach numbers and the low unit Reynolds number and 
fairly small size of most hypersonic facilities. The nozzle-wall boundary layer is, in 
most instances, a thick, readily accessible, turbulent shear flow. However, a number of 
problem areas arise when the genesis of this type of boundary layer is considered. The 
original turbulence structure is formed in a low-speed settling-chamber-wall flow. The 
boundary layer is next expanded to the throat region where high heat transfer is present 
and partial relaminarization often occurs. Farther down the nozzle, the main effect is a 
drop in external density until the final Mach number is reached. 

If a corresponding boundary layer had grown on a basic model in zero-pressure 
gradient, constant wall temperature, and at the same final Mach number, the flow would 
probably not be turbulent, at least for some of the available data (Reynolds number too 
low, see ref. 45). At low 6 + , probably the only reason the nozzle-wall boundary layers 
appear turbulent is the existence of a turbulent history in the low- speed settling- chamber 
region. Also, calculations carried out during the investigation reported in reference 1 
indicate that, from the throat to the exit of a Mach 19 hypersonic nozzle (assuming a 
fairly hot throat and a cooled test-section wall), the wall boundary layer thickens appre- 
ciably but does not entrain a significant amount of mass. Therefore, history effects in 
terms of distorted temperature profiles tend to remain with the flow rather than be 
"washed out" or relaxed. (See ref. 26.) 
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Therefore, on an actual vehicle where transition occurs naturally and in a high Mach 
number region, turbulent flow would probably only exist at fairly high 6 + values, and it 
would not be necessary to account for the increases in ( 4 ) and K which occur at 

low 5 + (figs, 12(a) and 12(b)), except perhaps in a method such as presented in refer- 
ences 5 and 11 which compute through transition. During transition, natural or tripped, 
a turbulence model should probably be corrected for low Reynolds number effects as 
well as intermittency. 

The ( — values for the outer region are shown in figure 12(c). These 

V u e 6 W 

\ e i/max , 7 \ 

results behave in a manner similar to the (-] values in figure 12(b). Most of the 


values at low 6 + are above the 0.016 to 0.018 eddy-viscosity range generally used in 
computational techniques. (See refs. 6, 7, 10, and 11.) 

As a possible alternate method of correlating the low- and high-speed "low Reynolds 
number effect," a Reynolds number was used based on sublayer thickness, wall density, 


u e y sl p w y sl u r w 

viscosity, and free-stream velocity; that is, R s j = — - , where -^±— 


10 was 


used to determine 


The K 


’ ft! 


, and 


values are shown in fig- 


max \ u e^i/max 

ure 13 as a function of R s j. Again, considerable scatter appears in the data, but the 
trends do not seem to be as well matched between the low- and high-speed data as in fig- 
ure 12 (except for K, fig. 13(a)); therefore, the use of 6 + as a correlating factor is 
favored. 


It should be noted that the nature of the variation of 


, u e^i/max 


, or K 


with 6 or R s | is probably a function of wall cooling. At the present time, insuffi- 
cient data exist to attempt to settle this question, and the possibility of wall-cooling 
effects is raised as a cautionary note in the use of the results shown in figure 12 or 13. 
More data at hypersonic conditions, especially with cold walls, are obviously needed 
before any definite conclusions can be made as to the correct scaling parameters for low 
Reynolds number effects at hypersonic cold-wall conditions. 


EVALUATION OF WALL-DAMPING EXPRESSION 


Hypersonic turbulent boundary layers generally have fairly thick sublayer regions 
because of the low density near the wall. In this near-wall region, turbulence damping 
occurs and the mixing-length distributions deviate from the usual l = Ky variation. The 
most commonly employed expression to account for this effect is equation (8), the Van 
Driest damping function (ref. 29). For compressible flows, there are several possible 
methods to use in evaluating r, /j. , and p appearing in equation (9), which is the 
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damping-scale expression. The most straightforward method is to evaluate all of these 
quantities at wall conditions. Based on the successful use of wall properties in refer- 
ence 41 to correlate the entire profile even where the density change across the boundary 
layer is a factor of several hundred, the use of wall conditions seems a reasonable 
approach. However, other suggestions, such as using sublayer averaged or local condi- 
tions (ref. 20) on p and p or using local shear (ref. 46), have been made. In the pres- 
ent case, the shear is approximately constant in the near-wall region (no large pressure 
gradient, wall blowing, or suction); therefore, the question of whether or not to use local 
shear cannot be investigated. However, especially for the Mach 20.8 investigation (fig. 2) 
where a sizable change in density exists across the sublayer, the question of whether to 
use wall or local p and p in equation (9) can be examined. Note that there are sev- 
eral other wall-damping expressions which could be used. (See refs. 7 and 16, for 
example.) 

Several of the profiles examined herein exhibit fairly extensive regions of wall- 
damping effects upon the l distribution. As an extreme case, the results shown in 
figure 6 indicate an extensive region near the wall with completely damped turbulence , 
a clue that the boundary layer may be "overdamped" as indicated in the previous section. 

By using the K values given in table I (obtained from figs. 2 to 11), several of the 
computed 1/5 distributions have been used in equation (8) to compute a variation of A + 
with y/8 for the cases where (1) p and p are evaluated at the wall temperature and 
(2) p and p are evaluated as functions of y. If the damping expression is correct 
and the 1/5 variation accurate, A + should be a constant with y/8. In addition, if the 
method of evaluation for p and p is correct, the value of A + should be close to 26, 
which is the incompressible level generally used in an impermeable-wall situation. (See 
refs. 1 and 29.) 

The resulting. A + variations are shown in figure 14, The values exhibit large 
deviations from the expected value of 26, probably due more to inaccuracies in the l/5 
variations reported herein than in the form of the damping function. By treating the 
deviation from the 26 value of A + as a random error, the root-mean-square deviation 
and average error were obtained for both methods of evaluating p and p in equa- 
tion (9). The resulting errors are shown in figure 15 as averages for all the cases of 
figure 14. From this figure, it is evident that (1) the value of 26 is somewhat too low, 
and (2) the method which provides the least error results from an evaluation of p and 
p at wall conditions. 


CONCLUDING REMARKS 

A method is developed to obtain shear-stress, eddy-viscosity, and mixing-length 
distributions from hypersonic turbulent boundary-layer profile data by using the 
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assumption of local similarity. Application of the method to published profile data indi- 
cates an increase in the maximum value of mixing length and eddy viscosity as the 
Reynolds number decreases. The level and trend of this low Reynolds number effect 
agrees with similar results previously obtained at low speeds when wall conditions are 
used in the correlating parameter. Examination of the wall-damping effect on the mixing 
length indicates more validity for evaluating the Van Driest damping scale by using wall 
conditions than by using local sublayer properties and that the impermeable -wall damping 
constant may be somewhat greater at hypersonic conditions than the usual low-speed value 
of 26. In several instances, areas of disagreement occurred between the expected and 
observed behavior of the data, usually concerning the shape of the published velocity pro- 
files. Indications are that the total-temperature data obtained in hypersonic, low Reynolds 
number, cold-wall boundary layers may not be suitably corrected for conditions near the 
wall. There is a clear need for more complete investigations of the mean flow structure 
of hypersonic turbulent boundary layers, especially at cold-wall conditions. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., June 30, 1971. 
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TABLE I.- EDGE AND PROFILE PARAMETERS 


Parameter 

Profile examined 

Axisymmetric 
nozzle wall 

Axisymmetric 
nozzle -wall 
extension 
(ref. 30) 

Conical 
nozzle wall 
(ref. 31) 

Two-dimensional 
nozzle wall 
(ref. 33) 

Hollow cylinder 

Flat plate 
(ref. 36) 

Ref. 21 

Ref. 32 

Ref. 34 

Ref. 35 

c i 

2 

9.75 X 10‘ 5 

7.50 X 10- 4 
4.00 X 10-4 

4.50 X 10‘ 4 

5.25x10-4 

2.60X10-4 

3.60 X 10-4 

4.3 X 10” 4 

3.00x10-4 

3.00 x lO' 4 

9.5 x lO -4 


6.90 x 10 3 

1.12 X 10 3 

8.00 x 10 3 

12.80 X 10 3 

57 X 10 3 

5.0 X 10 3 

1.4 X 10 4 

1.28 x 10 4 

2.53 x 10 3 

M e 

20.0 

12.00 

5.76 

10.65 

4.67 

4.7 

6.02 

5.98 

6.5 

( 6 \( dp ®l 

-7.0 X 10- 5 

-6.1 X 10-4 

2.2 X 10'4 

-4.7 X 10*4 

0 

0 

0 

0 

0 

(p e u e fteJ 

d6 

dx 

0.056 

0.038 

0.017 

0.024 

0.007 

0. 0085 

0.0116 

0.010 

0.030 

1m 

To 

1.000 

0.295 



0.700 


0.895 

0.450 


P W 

Pe 

0.007 

0.119 

0.1787 

m 

0.259 

0.275 

0.135 

0.274 

0.218 

dx 

0.06 

0.052 



0 

0 

0 

0 

0 

R sl 

84.73 

125.96 

172.48 



268.22 





5 

r w 

0.4735 

0.4670 

0.4400 

0.4930 

0 


1 


Jj 

5 + 

106.5 

66.5 

48.6 


.v.'M 

1 ■ 

2841.0 



1 


K 

0.65 




0.49 


0.61 

0.51 






/_L_\ 

0.060 

0.0132 


0.0162 



0,0172 

0.0135 

0.0320 

( u e 5 i) 

\ e 1 /max 



( L ) 

W max 

0.2 to 0.3 

0.075 


| 

0.068 to 0.078 

0.1 to 0.12 

0.06 to 0.08 

0.06 to 0.08 


5 ^ 

8.9945 

2.8109 

3.3905 

2.9219 

2.5624 

2.5484 

4.9275 


3.9572 


6 

0.0515 

0.1953 

0.1077 

0.1327 

0.1406 

0.1420 

0.1048 

0.0837 

0.1056 


















u 

u 
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Figure 1.- Velocity profiles from reference 21 at M e = 20 in assumed 

T 

similarity coordinates. ^ = 1.0; R 0 ~ 7 x 10 3 . 
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Shear stress 



layer at M e = 20.8. Profile data from reference 21; R fl = 6.90 x 10^ 











(b) Eddy-viscosity distribution 



Using computed shear 
distribution 



(c) Mixing-length distribution 
Figure 4.- Concluded. 





Figure 5.- Shear-stress, eddy-viscosity, and mixing-length distributions for a conical nozzle-wall boundary layer 

at M e = 10.65. Profile data from reference 31; R q = 12.80 X 10^. 



Present calculations 



Figure 5.- Continued. 





Figure 5.- Concluded. 






Figure 6.- Continued. 







boundary layer at M e = 4.67. Profile data from reference 33; R# = 57 x 10^ 
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igure 7.- Conclui 





(a) Shear-stress distribution. 

Figure 8.- Shear-stress, eddy-viscosity, and mixing-length distributions for a two-dimensional nozzle-wall 
boundary layer at M e = 4.7. Profile data from reference 33; Rg = 5.0 x 10^. 






Figure 9.- Shear-stress, eddy-viscosity, and mixing-length distributions for a boundary layer on a hollow cylinder 

at Me = 6.02. Profile data from reference 34; Rg = 1.4 x 10^. 








Figure 10.- Shear-stress, eddy-viscosity, and mixing-length distributions for a boundary layer on a hollow cylinder 

it M e = 5.98. Profile data from reference 35; Rg = 1.28 x 10^. 





Shear stress 



6.5. Profile data from reference 36; Rg = 2.53 X 10 
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(a) Mixing-length Prandtl wall slope. 

Figure 12.- Correlation of results from present turbulence model with low-speed data on a 6 + basis. 




(b) Mixing length in outer region. 
Figure 12.- Continued. 
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(c) Maximum eddy viscosity. 
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(b) Root-mean-square error. 
Figure 15.- Concluded. 
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